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The influence of nonperturbative fields on instantons in quantum chromodynamics 
is studied. Nonperturbative vacuum is described in terms of nonlocal gauge invari- 
ant vacuum averages of gluon field strength. Effective action for instanton is derived 
in bifocal approximation and it is demonstrated that stochastic background gluon 
fields are responsible for infra-red (IR) stabilization of instantons. Dependence of 
characteristic instanton size on gluon condensate and correlation length in nonper- 
turbative vacuum is found. Comparison of obtained instanton size distribution with 

O ■ lattice data is made. 

(N 

P^. PACS numbers: ll.15.Tk, 12.38.Lg 

o 



(N 
O 



(N 
> 

in 



o 

ix. 



INTRODUCTION 



Instanton is the first explicit example of nonperturbative (NP) quantum fluctuation of 
gluon fleld in QCD. It was introduced in 1975 by Polyakov and coauthors [|T|. An important 
^ ! development originated with 't Hooft's classic paper p|, in which he calculated the semi- 
Q \ classical tunneling rate. Instanton gas as a model of QCD vacuum was proposed in pioneer 
' works by Callan, Dashen and Gross |P, ^ . These topologically nontrivial fleld conflgurations 
give a possible explanation of some problems of quantum chromodynamics. Instantons allow 
to explain anomalous breaking of U{1)a symmetry and the rj' mass |^. Spontaneous chiral 
symmetry breaking (SCSB) can be explained with the help of instanton and anti-instanton 
fleld conflgurations in QCD vacuum [0]. An important role of instantons in scalar and 
^ . pseudoscalar channels was pointed out in paper [§]. Taking into account instantons is of 

crucial importance for many phenomena of QCD (see and references therein). 
k><( ; At the same time, there is a number of serious problems in instanton physics. The flrst 

\^ • is the IR divergence of integrals over instanton size p at large p. This makes it impossible 
^ to calculate instantons' contribution to some physical quantities, such as vacuum gluon 

condensate. Second, "area law" for Wilson loop can not be explained in instanton gas model, 
hence quasiclassical instanton anti-instanton vacuum lacks conflnement which is responsible 
for hadron spectra. It was recently demonstrated that Casimir scaling which is observed 
in lattice calculations for the interaction potential between heavy quarks QQ in different 
representations of 5't/(3) color group |10[ can be violated in instanton gas |]1T|. This matter 



is also discussed in |]T2[ 

There were many theoretical works aimed to solve the problem of instanton instability in 
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it's size p^. To some extent they all were based on an attempt to stabilize instanton ensemble 
due to effects of interaction between pseudoparticles |jl3|. Today the most popular is the 



model of " instanton liquid" , which was phenomenologically formulated by Shuryak ||T^ , 15 
Using hypothesis that instantons dominate in vacuum averages of local gludynamical opera- 
tors, and taking into account dipole-dipole interaction between instanton and antiinstanton, 
he determined main characteristics of the medium: average distance between pseudoparti- 
cles is ^ ~ 1 fm and their average size is p ~ 1/3 fm. Thus, a small quantity p/R ~ 1/3 
appears and a picture of "granular vacuum" emerges, i.e. vacuum consists of well separated 
and therefore not very much deformed instantons and anti-instantons. Nevertheless, action 
of pair interaction expl^j^f} ^ 1 appears to be important for the structure of instanton 
ensemble, and therefore one can make a conclusion that instanton ensemble in QCD is not 



a dilute gas, but an interacting liquid ||14 



Quantitatively similar results for parameters of instanton liquid were obtained by Di- 



akonov and Petrov [|16|. They have suggested that vacuum in gluodynamics consists of 
superposition of instantons and antiinstantons, and found that stabilization is related to 
discovered by them classical repulsion between pseudoparticles. However, further devel- 
opment [0 revealed that instanton ensemble can not be stabilized due to purely classical 
interaction. Thus, the mechanism for the suppression of large-size instantons in the ensemble 
of topologically non-trivial fields is still not understood. 

A more natural assumption is that there exist other nonperturbative fields apart from 
instantons in the vacuum, which allow to explain listed above problems. Interaction of 
small-size instanton p < 0.2 fm with long-wave gluon fluctuations, described by local vacuum 
condensate ((fi'G" )^), was studied in [0. It was demonstrated that those fields lead to even 



faster growth of instanton density with increasing of it's size p. On the other hand, in ||T9| , |20 



the study of instantons in stochastic QCD vacuum, parametrized by nonlocal gauge invariant 
vacuum averages of gluon field {tTG^u{x)^{x,y)G^\{y)(^{ii,x)) (here is a parallel 

transporter) was started. It was shown that standard perturbation theory undergoes a 
change in NP vacuum |T^, and that nonlocal interaction between large-size instanton 



(p > Ifm) and nonperturbative background stops the infrared inflation Next, it was 



demonstrated in that nonlocal interaction of instanton with NP vacuum fields stabilizes 
it on scales of order of correlation length in vacuum condensate. 

In the present paper we develop gauge invariant method for calculating effective action 
of an instanton in the NP vacuum. We show that instanton exists as a stable topologically 
nontrivial field configuration with characteristic size pc- Value of pc functionally depends 
on properties of bifocal correlator of NP field (tr G(x)<l'(a;, x)), i.e. it depends 

on two quantities: (G^) - gluon condensate, and "measure of vacuum inhomogeneity" 



9 



- correlation length in NP vacuum. We don't discuss the problem of instantons distri- 
bution in 4D-Euclidean space, i.e. the question about their density N/V. Consideration 
of this problem in the framework of our approach requires studying interacting instanton- 
antiinstanton ensemble against the background of nonperturbative fluctuations, described 
by gauge invariant vacuum averages of gluon field. 

This paper is organized as follows. Section H describes general formalism for effective 
action of an instanton in nonperturbative background. In Section |T| we study one-loop 
renormalization of instanton action in the presence of NP fields. It is shown that "perturba- 



Callan, Dashen and Gross suggested in their pioneer works y| to solve some of listed problems with the 
help of break-down of instantons into mcron-antimcron pairs. 



tive" part of effective instanton action is "freezed" in the NP background, and approaches 
constant value when instanton size is large enough p > ~ lGeV~^. Section |^ is 

devoted to the study of "direct" interaction between instanton and NP fields. NP fields 
are parametrized by gauge invariant nonlocal vacuum averages of gluon field (correlators). 
Using cluster expansion we find effective action for an instanton in bilocal approximation 
in a gauge invariant way. Results of numerical calculations and discussion of results are 
given in Section |V|. Derivation of one-loop renormalization of effective instanton action in 
nonperturbative vacuum in coordinate space, as well as some mathematical additions are 
brought out to Appendixes. 

II. GENERAL FORMALISM 

The influence of NP fluctuations on instanton can be separated into two parts. First, 
perturbation theory undergoes a change in NP background, which results in the change 
of standard one-loop renormalization of instanton action. Second, there appears a direct 
nonlocal interaction of instanton with NP background field. 

Standard Euclidean action of gluodynamics has the form 

S[A] = ^J d'xtrF'^M = \j d'xF^,AA]F^,u[A]. (1) 

where -F^i/[v4] = d^Ay—duA^—i[Afi, is the strength of gluon field and we use the Hermitian 
matrix form for gauge fields Af^{x) = goA1j^{x)t°- /2 and trt"t^ = 5"''/2. We decompose as 

A, = * + B, + a„ (2) 

where A^^^ is an instanton-like field configuration with a unit topological charge QtI^™^*] = 
1; ttfj, is quantum field (expansion in reduces to perturbation theory, which in gluody- 
namics leads to asymptotic freedom); is nonperturbative background field (with zero 
topological charge), which can be parametrized by gauge invariant nonlocal vacuum aver- 
ages of gluon field strength^. 

In general case effective action for instanton in NP vacuum takes the form 

Z = e-Seff [A'-] = j [^^^] ^g-5[A'-t+B+a] \^ ^ (3) 

where (...) implies averaging over background field B^, 

{0{B)) = I dp{B)d{B) (4) 

and dfi{B) is the measure of integration over NP fields, explicit form of which is not impor- 
tant for the following consideration. 



^ In operator product expansion method and in QCD sum rules nonperturbative field is characterized by a 
set of local gluon condensates (G^), (G^), . . . 



Expanding S[A] up to the second power in one obtains Z = {Zi{B)Z2{B)) , where 



[Da^] det(Vj)exp 



9o 



tr / d^x {-(V^a,)2 + 2iF^,[a^, a,]} 



Z2{B) = exp{-^[A"^^* + B] + 

Here we use the notation A = A™"** + B, F^^ = F^^[A], 
a covariant derivative. Integration over a and B in (|^),(| 



a/.,a^J| , (5) 
(6) 

and V^a^, = d^a^ - i[A^,a^] is 
]) corresponds to averaging over 
fields that are responsible for the physics at different scales. Integration over takes 
into account perturbative gluons and describes phenomena at small distances. Averaging 
over B^ (formally interaction with gluon condensate) accounts for phenomena at scales of 
confinement radius. Therefore it is physically clear that averaging factorizes and one obtains 

Z^{Z,{B)){Z,{B)). (7) 

It should be noted that in the limit of infinite number of colors oo factorization (J^) is 

exact Z{Nc oo) = {Zi{B)){Z2{B)). This allows us to write effective action of instanton 
in NP vacuum as a sum of two terms, "perturbative" and "nonperturbative" : 



Ses[A 



cNPr 



instl 



qNP 



ln(Zi(5)), 

- ln(Z2(5)) = - In (exp{-^[A^'^^* + B] + . 



(8) 
(9) 
(10) 



III. ONE-LOOP RENORMALIZATION OF INSTANTON ACTION IN NP 

VACUUM 

The general expression for one-instanton field configuration has the well known form 

- (^) . (11) 

where 

6=1,2,. ..iV^-l; /5 = 1,2,3. ^^^^ 

r/^j, - 't Hooft symbols, and matrix R'^^ ensures embedding of instanton into SU{Nc) group, 
R^^ satisfied relations 

In singular gauge profile function f{z) satisfies boundary conditions /(O) = 1, /(oo) = 
and the classical solution has the form 

The probability to find an instanton is determined by the classical action functional S'ci[v4] 
for the solution (|1|D; that is 

w ~ exp{-5ei[Af*]} = exp{-87rV^7o'} • (15) 



The preexponential factor was calculated in [0]. The result for one-instanton contribution 
to the QCD partition function is Zj = J dn{p,xo, R), where dn is the differential instanton 
density 



dn{p, xq, R) = [dR]d'^xo-jdo{p), (16) 
4.6exp{-1.68iVj / / \ 



7r2(iV,-l)!(iV,-2)! \g^ip)J ^ { g\p) 
In the two-loop approximation in gluodynamics, the coupling constant g'^{p) is given by 



6 In 



1 



hi 



In In 



1 



-0 



ln(l/pA) 
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iVc, &i = ?iVc, (18) 



A 



= Apy corresponds to the Pauli-Villars regularization scheme, A ~ 200 MeV. 

Callan, Dashen and Gross [HI were first to show that in a constant gauge field instanton 



behaves as a colored four- dimensional dipole. Shifman, Vainshtein and Zakharov gen- 
eralized this result to the case of interaction between a small-size instanton, p < 0.2 fm, and 
nonperturbative long-wave fluctuations that are described by the local vacuum condensate 
(G^). As a result, do{p) is replaced by des{p), where 



dM oc (Ap)" 1 + 



(iV2 - 1)^' 



(19) 



Thus, we arrive at the well-known problem of the infrared inflation of instantons in p. 

Let us now consider how perturbative renormalization of effective instanton action 
changes in NP vacuum |2^. We will make use of background field method, developed 
in [22, 23], and stated in a convenient way by Polyakov [23] when discussing asymptotic 



freedom in nonabelian gauge theories 
In the absence of the term 

_ for Zi. 

det-V^(V 



n m expression 



Fi_iu[a'iM, ttiy], there are four independent components of the 
Integration with respect to these independent compo- 

Taking into account the ghost determinant in (^, we obtain 



field a 
nents yields 

det"^/^(V2) Thus, there are all in all two physical polarizations for the field instead 

of four, as it must be. The last term in (|^), ~ Ff^i^[af^,ai^] describes the interaction of the 
external field with the spin of the gluon a^. 

Hence, there exist two effects. First, there is the motion of charged particles in the 
external gluomagnetic field F^^^. This motion leads to screening similar to orbital Landau 
diamagnetism. Second, there is the direct interaction of Ff^^, with the spin of the field a^. In 
the second order in the external field, this interaction yields an antiscreening term, which is 
analogous to that describing the Pauli paramagnetic effect. With logarithmic accuracy, these 
two phenomena can be separated. In the absence of the term ~ F^i^[a^,ai^] from (1), the 
expression for S"^ = — In Zi then completely corresponds to the effective action of massless 
scalar electrodynamics with allowance for the isotopic factor. Hence, in the second order in 
the background field, the " diamagnetic" contribution to can be written as 



cP 



(20) 



where Ildia is given by the diagrams in Fig. ^ 




Figure 1: Diagrams for Hdia 

Let us first consider a simpler case of A = ^4™^*; that is, we set the nonperturbative field 
-B^ = and demonstrate how the standard perturbative renormalization of the classical 
instanton action is reproduced. We then have 

K^M = 2r''f^U(^:\q) = 2NJ-'U^;:\q) (21) 

The factor of 2 in (^) stems from two physical polarizations of the field a^. The expres- 
sion for the polarization operator in massless scalar electrodynamics has the form 



KA^) = -4 y (2^ p^ip + qy + 2^- J (2^^ = n^^^ - ^22) 

The well-known result for Hi in the leading logarithmic approximation is given by 



1 . A, 



where Aq is the ultraviolet cutoff. Thus, we have 



nldia(?) = ^NJ'^'iqX, - q,q.)U,iq') (24) 



Considering that the linear component satisfies the relation 



KM)F%{-^) = - q,qu)Al{q)Al{-q) (25) 

and substituting (|^) and (pSj) into (|20|), we obtain 



sl. = \j ^ndia(g^)F;,(g)F;,(-g), 
nrf) = 4iv.n,(,^) = ^in^ 



(26) 



It should be noted that, in (^), the field strength involves only the linear component 
~ dfj_Ay. However, from gauge invariance and invariance under renormalization-group trans- 
formations, it follows that cubic (~ A'^) and fourth-order (~ A^) terms appear in such a 
way that they supplement F^i^ to the standard non-Abelian form. 



Following Polyakov [Q, we represent the paramagnetic term in the form 



Ag^J {2ny ^' J {2^^ p{p + q) 

1 f dS 



J ^np_(g^)F;,(g)F;,(-g), 



where 



iv. . a: 



2 



np..a(g^) = -^in-f (28) 

Bringing the above results together, we find that one-loop renormalized instanton action is 
given by 



^ - 4 y (2^ U ^ "^'^ ' ^;.(^)^;.(-^) (29) 



where 

n(g2) = ndia(g') + npara(g') 

= ^fl_lVn^ = -iliV^ln^ (30) 

As might have been expected, involves the effective charge g'^{q) defined at the external 
field momentum q 

g^{q) gl 3 q^ ^ ' 

Using the standard renormalization-group arguments, we perform normalization at the mo- 
mentum 

We then have 

We consider assuming that the background field is an instanton. Because is a gauge- 
invariant quantity, we use the regular gauge 

i^;.(^) = -4^^. ^^^^ 

Considering that the Fourier transform in four-dimensional Euclidean space is given by 

r 47r2 
f{q) = / c/W«V(a;) = / r^drf{r)J^{qr), r = \x\ (35) 



Q Jo 



we obtain 



(36) 



'qJo (r2 + p2)2 



where Kq{x) is a McDonald function. It can easily be shown that the classical action for 
solution (^6D has the standard instanton value 



rq'dqKl(pq)^%. (37) 

9q Jo 9o 

Substituting (^) into (^9[) , we arrive at 

= STT^Sp^ / q^dqK'.ipq)^ (38) 
Let us introduce the dimensionless variable x = pq and rewrite (|38D in the form 

^PMinsti 8^23 / AXg(a;) , (39) 

7o r(a;/p) 

where 

1 6 1 ^ 1 2 

^2(x/p) " 16^ AV " 16^ AV ^ 16^ """^ 
Thus, we finally obtain 

S^(p) = ^, , + const = 61n h const . (41) 

9^{p) Ap 

We have reproduced the standard expression for S^{p). In the instanton density d{p), the 
constant from (^) appears as a preexponential factor, which is immaterial in our approach. 
We will now study the modifications that arise in the above equations as the result of taking 
nonperturbative vacuum fluctuations into account. 

It was shown in that, for the polarization operator in the nonperturbative vacuum, 
the exponential decay of the form ~ |x — y\~^ exp{—m^\x — y\} is obtained in the infrared 
region instead of the standard power behavior ^ \x — y\~^ . Here ^m^, ~ 0.75mo++ ~ IGeV, 
and 0"^"^ is the lightest glueball (see Appendix We will use this fact in analyzing instanton 
in the nonperturbative vacuum. It can be seen from (|^) that averaging is performed over 
the field B. This means that the polarization operator in (^) is represented by the diagrams 
in Fig. |. 

Hence, in the expression ( p9D for the effective action, we must replace the perturbative 
expression for the polarization operator by the polarization operator Hnp(5'^) in the nonper- 
turbative background field. The latter was calculated in |24| and is given by 

n-('^)-T^=I^'"¥ (42) 

Effective action (^) is then written as 



^inst 
• 





^inst 
• 




^inst ^inst 
« » 




Figure 2: Diagrams for polarization operator IInp in the nonperturbative background field 



where 



Substituting the Fourier transform ( ^6|) of the instanton field into ( ^3|) and introducing the 
variable x = pq, we obtain 



247r' 



X ' doc -^^Q ( ^ } 2 



1 



b , + mlp^ 



(45) 



167r2 



In 



Because In is a slowly varying function and because the main contribution to the integral 
in (|45| ) comes from the region x ~ 1 {Kq{x 3> 1) — \Jt: /2xe~^), we can take l/5f^p(x) out 
of the integral with logarithmic accuracy and set x = 1. In this way, we obtain |jl9|, ^ 



h l/p2 + ml 



(46) 



From (^61), it follows that the perturbative expression is recovered for small-size instantons 
(p <^ 1/m*) and that, for p > 1/m*, we have 5^ const. 

Derivation of formula (|46|) in coordinate space is given in Appendix 0. 



IV. INTERACTION BETWEEN INSTANTON AND NONPERTURBATIVE 

VACUUM GLUON FIELDS 



We consider effect of NP fields on instanton, i.e. we evaluate (^'2). In this work we make 



use of the method of vacuum correlators, introduced in works of Dosch and Simonov [|25 



NP vacuum of QCD is described in terms of gauge invariant vacuum averages of gluon fields 
(correlators) 

= (tr (xi ) $ (xi , 0:2 ) ,.2 (x2 ) . . . Gf,^ (x„ ) $ (x„ , xi ) ) , 

where G^^ is gluon field strength, and $(x, y) = Pexp(i B^dz^) is a parallel transporter, 
which ensures gauge invariance. In many cases bilocal approximation appears to be suf- 
ficient for qualitative description of various physical phenomena in QCD. Moreover, there 
are indications that corrections due to higher correlators are small and amount to several 
percent ^ . 



The most general form of bilocal correlator, which follows from antisymmetry in Lorentz 
indices, is given by 



^^^^ {Sf,p6,ya - S^aSup) + ^^^{nf,np6ua + n^n^S^p - n^naS^p - n^Updpa) \ , (47) 



where 



Gp„{x, xo) = $(xo, x)Gpu{x)^{x, xq), (48) 

^p, = Zp/\z\ = {x — y)p/\x — y\ is the unit vector, (G^) = {g^G'^^^^G'^^J) and, as it follows from 
normalization, D(0) + D(0) = 1. 

Connection between functions D{z), D{z) and standard functions D{z), Di{z), which are 
used in works ^ , is discussed in Appendix 0. 

Functions D{z) and D{z) have both perturbative and nonperturbative contributions. We 
will consider only nonperturbative part, because perturbative one has already been taken 
into account in S^g. Most data about nonperturbative components of these functions come 
from numerical simulations on lattice. Gluon condensate (G^) is also determined from 
lattice calculations, but there exists a widely used estimate for this value based on char- 
monium spectrum analysis and QCD sum rules ~ 0.5 GeV^ |9[. According to the 

lattice data D{z) and D{z) are exponentially decreasing functions D{z) = exp(— z/T^), 
D{z) = Aizexp{—z/Tg)/Tg, where Tg is the gluonic correlation length, which was measured 
on the lattice |2^, ^ and estimated analytically |3ll] to be Tg ~ 0.2 fm. Besides, according 



to lattice measurements Ai <^ Aq {Ai ~ Aq/10). Lattice data from paper of Di Giacomo ||2l 
are presented in table |I[ SU{3) full stands for chromodynamics with 4 quarks, while SU{2) 
and SU{3) quenched stand for pure SU{2) and SU{3) gluodynamics, respectively. Bilocal 



correlator was also measured on the lattice without using cooling method and correla- 
tion length was found to be of around 0.1 — 0.2 fm in quenched approximation. 

Gluon condensate (G^) = 1.02 ± O.lGeV^, obtained by Narison ||3^ basing on sum rules. 



is in good agreement with lattice calculations [^, {G^) = 0.87 GeV*, for SU{3) QCD. 

Table I: Lattice data [^] for correlation length 



Tg, fm 



SU{2) quenched 0.16 ± 0.02 
SU{3) quenched 0.22 ± 0.02 
SU{3) full 0.34 



To evaluate S^q defined in ( p!0| ) we use the cluster expansion, which is well known in 
statistical physics PB| : 

{exp{x))=exp(^^y (49) 



where (x) = ((x)); (x^) = {{x^)) + {x)^; (x^) = ((x^)) + 3{x){{x'^)) + (x)^- 



One can easily find coefficients in cluster expansion, rewriting it in the following way 



(1) = ((1)) 
(12) = ((12)) + (1)(2) 

(12 3) = ((12 3)) + (1)((2 3)) + (2)((13)) + (3)((12)) + (1)(2)(3) ^ > 



Next, we modify expression (p!0| ) for S^^ by adding S[B] to the exponent: 

'5eT[^^'''*] = - In (expl-^iA'-^^^ + B] + + S[B]}) . (51) 

In bilocal approximation this is equivalent to a change of normalization of partition 
function,^ which is not important for the following consideration. 

Nonperturbative part of the effective instanton action S'^q takes the form 

S^i = {S[A'''''' + B]- S[B] - S[A'°^*])- 

We don't go further than | (( {S[A'''^^ + B]-S[B]- ^[A'"^*^']) ^)) term. Next terms of cluster 
expansion decrease as 1/g'^"'. It is shown in that g'^{pc) ~ 4 on the scale of instanton 



typical size, and therefore 1/(7^ is a small parameter. Another small parameter is l/N^. 

By using Fock-Schwinger gauge x^A^ = x^B^ = (instanton field A^^^*' satisfies it owing 
to the properties of 't Hooft symbols) and taking into account that 

^[^inst ^ B]-S[B] - S[A'''''] = ^J rf'a;tr|-([Af B,] - K'^^S B,]f 

+ 2F^4A'-'']G^4B]-At{F^4A'-'']+G,,AB])[Af'\B,]] (53) 

one gets 

~ ^dia. + + 'S'para + "S*! + 5*2, (54) 

where in bilocal approximation 

^dia = / [([^r*' Bu] - [A'r, B,]Y] ) , (55) 

5para = / d'xd'y {ti {F'^\x)G ,.{x)) tr {F';:\y)G,M)) , (56) 

51 = j,j d'xdS (tr (i^r [Af S 5.])^ tr [Af \ B^])^) , (57) 

52 = ^ 1 d'xd'y (tr {Fp'G,^)^tT [Af \ B^])^) . (58) 



^ Indeed, starting from ( p^ ) one finds that in bilocal approximation (|5^) gets additional terms of the type 
(G^) / d'^x, which do not depend on ^4™''* and therefore result only in renormalization of Z2. 



We use notations ^dia (diamagnetic) and Spara (paramagnetic) for contributions 



and (|56[) into interaction of instanton with background field. Physical motivation for this is 



discussed in detail in 20 



We also take into account that 

{ii{G,MTBu - 5.4°^*))) = i^-^r''{Gl,Bl)A';:^'' cx rH'-' = 0. (59) 
In Fock-Schwinger gauge 

1 

B„W^../a.a.G„„(..). (60) 



and thus we substitute vacuum averages by correlators in (|55D-(p8D. In general form instan- 
ton field has the form 

*(x) = ^x,Xo)Af^ix - Xo)$(xo,x). (61) 
/ - \ 

Here $(a;, y) = Pexp i J B^dz^ j is the parallel transporter, and A^™^{x — Xq) is instanton 



in fixed (singular) gauge (|TT]),(14). We then obtain for field strength 



= d,^r - 5.4°^' - ^u""'] = H^, ^o)i^;r(a: - Xo)<f (xo, x) + ... (62) 

Dots denote terms which contain background field B, and therefore lead to higher correla- 
tors and we take into account bilocal correlator only. 

Parallel transporters along with background field form gauge invariant combinations, 
which are substituted by correlators after averaging. For instance, one gets for iSpara 

^para = -7^ / d^xd% {ti {F^,{x)G ^,{x)) ti {F,^{y)G ,^{y))) = 



2g^ 
1 



2^4 
1 



2^7 



d'^xd^i tr ($(a;, Xq)F^^^{x - Xo)$(xo, x)G^^{x)) x 

X ii{^{y,xo)Ff:^{y-x,)^{xo.y)G,„{y))) = (63) 
d'^xd'y (tr {F^;^^ix - xo)G,.(x, xo)) tr {F^^^iy - x,)G,^{y, y,)) ) = 



^ f d'xdS {F^:^{x - x,)Y {Ff^^iy - xo))\g^Gl,{x, x,)G\Ay: ^o)> 



Analogous calculations for other contributions to S^^ (ill), lead to the following expressions 
(for the sake of simplicity we omit further index 'sing' for instanton field in singular gauge, 
AT^ = A.,). 



1 



^dia=- / d^x / ada / f3d(3r''^f^^{A;^Af^x,x,{g'Gliax,xo)G;,{(3x,xo)) 





-A^AtxxXpig'Gliax, xo)G;^{Px, xq))) 



1 



para 



8g' 



d'xd'yF",, {x)Fl {y) {g^G^, (x, x,)G' (y, xo) ) 



^1 



d'xd%Fl{x)A'{x)F';MAl{y)r'''f'^x 



1 1 



X / ada I pdPx^yr^{g'^G'i^{ax,Xo)GUPy,Xo)) 



(64) 



(65) 



(66) 







^^ = -Ygj '^''^'^'y J «^«C(^)^pa(?/)/'''^p(2/)l/€(^?'G;:,(x,xo)Gt(«^'^o)) (67) 



using representation (^) for gauge invariant condensate ((^^G" (x, xo)G^ {y, xq)) one gets 



dia 



12 m-1 







S, 



para 



48^2 _ 1 



d'^xd'^y 



F;^{x)D{x ~ y)Fl{y) + 



flU \ 



+ 4^^^^M^^F^(x)D(x - y) F^(y) 



(x - yf 



(69) 



24 ^"2 _ 1 



1 1 



abe rede 







(70) 



X FUx)Ai{x)F:MA%y)D{ax - Py) + 0{D) 



{G') 1 



d^xd%r"'FUy)Aliy)F:,ix)y, I adaD{x - ay) + 0{D) (71) 



I2g - 1 



We have taken into account, that 

x^Bf,{x) = Xf, ($(xo, x)5^(x)$(xo, x)) = $(xo, x) {x^Bf,{x)) <l>(xo, x) = 0, 



and thus 



(72) 
(73) 



B^{x) = Xjy y adaG f^iy{ax) ~ J 0!daG^„{ax, xq) + . . . , 
where in analogy with (|62D 
Gfiu{x) = d„B^ - d^B^ - i[B^, B^] = $(xo, x)Gf,^{x)^{x, Xq) + . . . = G^^{x, xq) + . . . (74) 



Thus we have obtained effective action for instanton in NP vacuum in bilocal approx- 
imation^. It can be seen from numerical calculation that typical instanton size in QCD 
is Pc ~ 0.25 -h 0.3 fm. For such p instanton field in it's center is strong F^^{x = Xq,Pc) = 
192/ p"^ ^ (G*^), and therefore classical instanton solution is not much deformed in the region 
\x\ < p, which gives the main contribution to integrals (|68|)-([7T|). 

The problem of asymptotic behavior of instanton solution far from the center ^ p 



was studied in detail in refs. |]T6|, B3, |3H|. Our numerical analysis shows that the value of 



Pc is almost not affected by the asymptotic of classical instanton solution provided that 
condensate (G^) and correlation length Tg have reasonable values. 

Next we have to define function D{z). We use gaussian form D{z) = exp{—p'^z^), 
p = l/Tg. D{z) is a monotonously decreasing function with typical correlation length 
Tg. Numerical calculations show that explicit form of D{z), as well as taking into account 
D{z) almost do not affect the value of instanton size in NP vacuum pc (see Appendix P). 

We use the standard profile function / = p^ /{p^ Certainly, exact instanton profile 

(i.e. profile that minimizes action) is different from this one. Nevertheless, knowing the 
dependence of S'eg on p we will be able to make a conclusion about instanton behavior 
(inflation or stabilization) and to determine it's typical size in NP vacuum. 

Thus we obtain after integration over space directions 

oo 

^'- = ^1^N[^1^ y^^(^^^V^(^) (75) 





oo 



_ -IGvr^ {G') 1 f xVe-(^'+y'U,{2xy) 

-^Para " ^4 ^4 ^2 _ ^ ^ ^ ^''^^V (^.2 + ^2)2(^2 + ^2)2 ^"') 



00 

^ -647r4 1 ,8 /" , , 



1 1 

X [ da [ dl3e-^''''''+^'y'\h{2al3xy) + h{2al3xy)) 



(77) 







,2 





where /„(x) = Jn{ix) is Bessel function of imaginary argument (Infeld function), ( = pp 
and 

1 1 



X 

where $(a;) = / is error function. 





^ The tensor structure of instanton was used to derive (|68|)-([7l|). 



We studied asymptotic behavior of instanton effective action and it's dependence on 
dimensionless parameter ( = ftp. For a small-size instanton p <^ 1/ p {( ^ 0) we found 
for 

5dia - -^2j;^(G')pUn{pp), (80) 

5para + S^ + S2^ -Kw-^P' + O {{G')p\pp)) . (81) 

In the opposite limit of large p ^ 1/p ^ oo) one has 

71-7/2 p3 

^dia--^^^(G^)-, (82) 

•S'para + Si + S2 ^ COUSt . (83) 

Differential instanton density is proportional to dn/d'^zdp oc exp(— S'efj). Thus, "diamag- 
netic" interaction of instanton with NP fields leads to its infrared stabilization in size p. 



V. DISCUSSION AND CONCLUSIONS 

In this work we have studied instanton behavior in NP vacuum, which is parametrized by 
gauge invariant vacuum averages of gluon field strength. We have derived effective instanton 
action in bifocal approximation and have shown that "diamagnetic" term S'dia leads to IR 
stabilization of instanton. 

Numerical results for Scs are plotted in Fig. |^. Three lines correspond to S'^(p), S^^{p) 
and S'efr(p) = S^^^p) + S^^{p) in the case of SU{?>) gluodynamics, at (G^) = 1.0 GeV^ and 
Tg = 0.3 fm. 

In Table |I| values of pc for different {G^) and Tg in case of SU (3) gluodynamics are listed. 
Values of pc in case of QCD with Nf = 2 are presented in Table [IT^ 

Instanton differential density dn/d'^zdp ~ exp(— S'efr) and corresponding lattice data ||36 



are presented in Fig. ^ All graphs are normalized to the commonly accepted instanton 
density lfm~^. Instanton size p ~ 1/3 fm was first obtained in numerical calculations 



on the lattice in It should be mentioned that over the last years the lattice results 



have not much converged (compare the conferences @). Different groups roughly agree 
on instantons size within a factor of two, e.g. p = 0.3, . . . , 0.6 fm for SU{3) gluodynamics. 
There is no agreement at all concerning the density N/V. As a tendency, lattice studies 
give higher density and larger instantons than phenomenologically assumed. 

Thus, our results for pc are consistent with phenomenological value p and lattice data. 
Furthermore, we can present physical arguments to explain existing deviations with some 
lattice groups, which give larger than phenomenological values of instanton size. Lattice 
calculations include cooling procedure, during which some lattice configurations of gluon 
field are discarded. This procedure can result in a change in gluon condensate (G^), and thus 
instanton size distribution is calculated on the lattice at a value of gluon condensate (6*^)0001 
which is different from physical value (G^). It should be noted that uncooled instanton size 
distribution was studied in the work where authors introduced a new scaling variable 
("cooling radius"), which helps to extract information on the uncooled distribution. 



We show dependence of pc on (G^) for several values of Tg in Fig. ||. One can see that 
increase of (G^) results in decrease of instanton size, and that effect is a result of nonlocal 
"diamagnetic" interaction of instanton with NP fields. Fig. ^ shows function of Tg 

for several values of (G^). The smaller Tg, the larger instanton is. It is physically clear that 
less correlated NP fields {Tg 0) have smaller influence on instanton field configuration, 
which occupies 4D-Euclidean volume with characteristic size p {p ^ Tg). On the other hand, 
perturbative quantum fluctuations tend to inflate instanton, and therefore pc increases with 
the decrease of Tg. 

We did not go beyond bilocal approximation in this work. As mentioned above, this 
approximation is good enough not only for qualitative, but also for quantitative description 
of some phenomena in nonperturbative QCD [^. In the problem under consideration there 
are two small parameters. These are l/g'^{pc) ~ 0.15, . . . , 0.25 and 1/Nc, and there powers 
grow in each term of cluster expansion. Moreover, one can make an estimate of the leading 
terms in cluster expansion. It follows from our numerical calculations that in the region 
p ~ Pc dominating terms are S'dia, ^S^ia: l^dia. ^o on. Summing them up, one gets 
'S'dia + |5'^ia + ^^dia. + ••• = — ln(l — S'dia)- Heucc, taking into account next terms of cluster 
expansion not only conserves IR stabilization, but even leads to some decrease of pc. So, 
proposed model describes physics of single instanton stabilization in NP vacuum, not only 
quahtatively, but also quantitatively with rather good accuracy. 

Acknowledgements 

We are especially grateful to Yu.A. Simonov for stimulating discussions and useful re- 
marks, and we thank E.-M. Ilgenfritz, A.B. Kaidalov, V.A. Novikov and A.V. Yung for 
discussions of results. The financial support of RFFI grant 00-02-17836 and INTAS grant 
N 110 is gratefully acknowledged. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 

p, fm 

Figure 3: Effective action 5^ (dotted line), S^^ (dashed line) and Seg = 5^ + (solid line) 
{Nc = 3, Nf = 0, Tg = 0.3 fm, (G^) = 1.0 GeV^.) 
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Figure 4: Instanton density dn/d'^zdp and lattice data [^] {Nc = 3, Nj = 0, Tg = 0.2 fm, 
= 0.5 GeVl) 



0.35 



. 3 



0.25 



0.2 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 

GeV^ 

Figure 5: Instanton size as a function of gluon condensate {Nc = 3, Nf = 0) at Tg = 0.1 fm (dashed 
line), Tg = 0.2 fm (solid line), Tg = 0.3 fm (dotted line) 
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Figure 6: Instanton size as a function of correlation radius {Nc = 3, Nj = 0) at (G^) = 0.5 GeV^ 
(solid line), (G^) = 0.87 GeV^ (dotted line), (G^) = 1.0 GeV^ (dashed line) 



Table II: pc (fm) for different values of Tg and (G^), Nc = 3, Nf = 



GcV^ 


0.5 1.0 1.5 2.0 


Tg, fm 




0.1 


0.305 0.256 0.230 0.213 


0.15 


0.282 0.236 0.213 0.197 


0.2 


0.268 0.226 0.204 0.189 


0.25 


0.256 0.219 0.198 0.184 


0.3 


0.254 0.215 0.194 0.181 


0.35 


0.250 0.212 0.191 0.178 



Table III: pc (fm) for different values of Tg and (G^), = 3, Nf = 2 



{G^), GeV^ 


0.5 1.0 1.5 2.0 


Tg, fm 




0.1 


0.308 0.258 0.232 0.215 


0.15 


0.284 0.238 0.215 0.199 


0.2 


0.270 0.227 0.205 0.191 


0.25 


0.262 0.221 0.199 0.185 


0.3 


0.256 0.216 0.195 0.182 


0.35 


0.252 0.213 0.192 0.179 



Appendix A: DERIVATION OF S^^ IN COORDINATE SPACE 

Let us consider derivation of (|46|) in coordinate space |19|. 



We start with expression for partition function, and normalizing it to the case A™*^* = 
we obtain effective instanton action in nonperturbative background 




s!AAn+s^nAn (m) 

CP c 1 /'detV2 /detZ}2y/2\ 

'5eff = 5ci-ln — XTT^ (A2) 



where 



det VdetDg 

S^i = - In (exp{-5[A'°^* + B] + S[A''''''] + S[B]}) (A3) 
S'ci = 8n'^/gQ, 

= (y^[A"''' + B]f, (A4) 
(D')^, = + 2tF,4A'-'' + B] 

while subscript '0' implies that one should put A™*^* = 0. 

Keeping in (|A2D terms quadratic in instanton field and using (|6ll) one has e.g. for the 
ghost contribution (gluon contribution can be obtained in the same way, final result will 
contain, as always, two physical gluon polarizations) 

S^^ = J d'xd^ytTAf^ix)d,dM^,y)A:^^iy) (A5) 

where explicitly gauge-invariant two-point function 11 is given by 

Il{x,y) = ($adj(0,a;)Vo'(x,2/)<l>adj(l/,0)Vo-'(2/,a;)), (A6) 

here Vq^(x,?/) = {x\l/Vl\y) is the ghost propagator. One can distinguish in ( [AG] ) a two- 
gluon correlator with an additional condition that adjoint parallel transporter $adj is passing 
through the instanton position Xq = 0. The mass spectrum of IT is bounded from below by 
the two-gluon correlator without that condition IIo, and the latter has as the lowest state 
0"'"+ glueball with the mass around 1.5 GeV. In the large- A'^c approximation Ilo{Q) has only 



poles p3 



oo ^ 

Mq') = Yl 2 K'=47raadjn + const (A7) 

n=l " 



as in qq case. The high-excited states in ( [A7| ) yield 



no(g') ~ ^ ( + const, ^(x^ - ^'^"^ 



ml J ' ' ' r(x) (A8) 

mg = 47r(Tadj^M2(0++) 



Replacing 11 by IIo in (^1) ^ind performing Fourier transform one has 



where we use the notation F^,^ for instanton field 



One can approximate Ilo{q^) in the whole Euclidean region by |T^, |20|, |24 

no(g^: 



487r2 ^ A2 ' 



where is connected to rriQ through a relation 



2 2 —C 



0,56ml, C 



-^d) = 0,577. 



(AlO) 



(All) 



Inserting in {\K^) F":^^{q) = — Svr^r/" p^fro(pg), where Ko{z) is the McDonald function yields 



5P(p) = 247rM x'dxK'.ix) — 
Jo 9 { 

where we have defined the effective charge 



x) g^iX) 



(A12) 



1 h x'^ + mlp^ 11 

■In ——: — , 0= — Mr 



g'^{x) IQtc^ K^p^ 



(A13) 



which is obviously the same result as (^q) . 



Appendix B: BILOCAL CORRELATOR 



The commonly used parametrization of bilocal correlator 26, ^ 28 1 is different 
from (147|). It has the form 



dDi 



~^ i^Zft^p^va ~l~ ZpZ(j5 ZpZijSi/p ZijZpS p(j) 



dz' 



Clearly, one can establish correspondence between our definition and (|Bl| 



(Bl) 



D{z) + D^{z) ^ D{z) 
1 ^dDi — , 



2^ dz^ 



D(z) 



(B2) 



Correlator (|B1| ) was measured on the lattice ||2^. Functions D{z) and Di{z) were approxi- 
mated as 

D{z) = Aoexp (-l^l/T,) + ^ exp (-|2;|/A) 

(B3) 

Diiz) = A, exp (-l^l/Tg) + 4 exp (-|^|/A) 



Second term in the right hand side of equations ( p3|) corresponds to perturbative contribu- 
tion. Numerical values for Tg, obtained in are presented in the Table |. Moreover, it 
was shown in these works that ratio Ai/Aq is a small quantity, Ai/Aq ^ 0.1. Turning back 
to our notation {D and D), we get for nonperturbative part 



D{z) = iAo + A,)e-'/^^ 

D(Z): 



(B4) 



4T„ 



It is clear that in the region z Tg, which gives the main contribution to all quantities, 
calculated in Section D{z) <C D{z) because Ai -C Aq. 
In our calculations we use gaussian parametrization 



A,e-^'/^^ 



D{z) = Aoe 
DAz) = A,e-'"/'^" 



(B5) 



Thus, we have 



D{z) = iAo + A,)e-^'/^' 
1 z^ 



D(z) 



A^e-^'in 



(B6) 



and still D{z^ <^ D{z) in the region z ~ Tg. 



Appendix C: DEPENDENCE OF 5dia ON FUNCTIONS D AND D 



Let us consider numerical changes in effective action due to taking into account function 
D{z) in bifocal correlator (^71) . As stated above, nonlocal "diamagnetic" interaction leads 
to instanton stabilization in size. Correspondingly, we will consider changes of S'dia due to 
D and to talcing exponential parametrization of bifocal correlator. Starting from (|68|)-(|7TD 
and taking into account (|B^ ) one gets for S'dia (instead of ([75D): 



dia 



X 



/i4 N^-r J (a;2 + C2)2 



} } , . (CI) 

<f{x) = ada I3dl3 [(Aq + Ai)e-("~^)'"' - Ai{a - /?) Ve-^""^)'"' 





Aq ~ 0.9, A 



0.1 



Fig. 1^ shows S'dia calculated according to formulas ([f5|)-([79|) and (|CT|). One can see, 
that the difference is very small, and therefore our approximation, when we neglect D{z) 
as compared to D{z) is justified. In the same figure there also is shown graph for Sdia, 
calculated with D{z) = e^'^'/^s, D(z) = 0. It is clear that exact form of function D{z) does 
not significantly influence the result. The only property this function should have is to be 
monotone decreasing function of z with characteristic correlation length Tg. 
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Figure 7: 5dia, with (dashed line) and without (solid line) taking D{z) into account. Dash-dotted 
Une - S'dia at D{z) = exp{-|2;|/rg}. 
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